arXiv: 1509.00927vl [physics.plasm-ph] 3 Sep 2015 


Modeling of Inelastic Collisions in a Multifluid Plasma: 

Excitation and Deexcitation®^) 

Hai P. and Jean-Luc Cambier^' 

Department of Mathematics, University of California, Los Angeles, 

California 90095 

Air Force Research Laboratory, Edwards AFB, California 93524 
(Dated: September 4, 2015) 

We describe here a model for inelastic collisions for electronic excitation and deexci¬ 
tation processes in a general, multifluid plasma. The model is derived from kinetic 
theory, and applicable to any mixture and mass ratio. The principle of detailed bal¬ 
ance is strictly enforced, and the model is consistent with all asymptotic limits. The 
results are verified with direct Monte Carlo calculations, and various numerical tests 
are conducted for the case of an electron-hydrogen two-fluid system, using a generic, 
semi-classical model of collision cross sections. We find that in some cases, the contri¬ 
bution of inelastic collisions to the momentum and thermal resistance coefficients is 
not negligible, in contrast to the assumptions of current multifluid models. This fun¬ 
damental model is also applied to ionization and recombination processes, the studies 
on which are currently underway. 

Keywords: Non-equilibrium plasma, collision, multifluid, Collisional-Radiative model 


Distribution A. Approved for public release; distribution unlimited 
Corresponding author; Electronic mail: hai.le@ucla.edu 
‘^^Electronic mail: jean_luc.cambier@us.af.mil 


1 



I. Introduction 


Modeling of nonequilibrium processes in a low-temperature partially ionized plasma is of 
particular interest to a wide range of technical fields such as gas discharge, electric propulsion, 
spectroscopic and laser diagnostics, and material science^^^. The complexity of the model is 
largely due to the characterization of various collisional and radiative processes, occurring at 
a wide range of spatial and temporal scales^’^. Although the fundamental physical processes 
may be individually known, it is not always clear how their combination affects the overall 
operation, or at what level of detail this process needs to be modeled. The current state of 
the art for modeling detailed chemical kinetics of a low temperature plasma is the collisional- 
radiative (CR) model, hrst proposed by Bates et ah in 1962®’^. CR models are now commonly 
used in studies of plasma discharge, plasma-assisted combustion, and hypersonics®'^^. The 
advantage of a CR model is two-fold. First, strong deviations from equilibrium of the internal 
states can be captured accurately when CR models are employed. In addition, ab initio 
cross section data can be directly incorporated in the CR model, leading to a very accurate 
prediction of the thermochemical kinetics of the system. 

There are, however, two issues with this modeling approach. The hrst arises from the 
complexity of the physical processes needed to be captured in the model. The required 
level of detail of the CR model is typically not known a priori and is possibly changing in 
a dynamical fashion as the system evolves in time. This can be resolved by coarse-graining 
techniques, which reduce the complexity of the kinetics to avoid solving a large system of 
equations^^^^®. The second issue comes from translational nonequilibrium, often found in a 
discharge, where we have both a bulk plasma (continuum) and a highly energetic component 
(kinetic), e.g., electrons emitted from the cathode. A proper treatment of this energetic 
beam-like component requires extending the solution of the CR kinetics to the so-called 
non-Maxwelllian regime^^^^®. These simulations are typically very expensive and therefore 
limited to zero- or one-dimensional systems. In addition, space charge effects within the bulk 
plasma can also become important, requiring further separation, i.e., ions and electrons. A 
natural solution to this problem is to use a hybrid approach, decomposing the system into 
a continuum and a kinetic component^^. The continuum component can be solved by fluid 
equations, and the kinetic component by (for example) a particle method. Although the 
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idea seems quite intuitive, proper treatment of the coupling between the fluid and kinetic 
components is highly non-trivial. There appears to be no unique and consistent coupling 
methodology, and the choice is highly problem-dependent. This coupling issue is currently 
being addressed by many researchers and is outside of the scope of the current work. 

We are concerned here with an alternative approach, the so-called multifluid model, which 
decomposes the plasma into several fluid components. For example, in a discharge configu¬ 
ration, one could have 4 different fluids, namely the neutrals, ions, bulk electrons, and the 
energetic electrons. The only required assumption is that collisions among particles within 
the same fluid are sufficiently fast to maintain a Maxwellian distribution. The validity of 
this assumption is not always well known. Nevertheless this approach is attractive, since 
it is much faster than a fully kinetic solver, and unambiguous since at the same time, the 
approach relies on kinetic theory for the treatment of coupling terms between different ffuids. 
Furthermore, one can rely on fast implicit methods to solve these fluid equations and exam¬ 
ine long time behavior of the system; this offers a big advantage over kinetic simulations, 
which are only suitable for problems with short time scales. Multifluid models are commonly 
used in simulations of astrophysical plasmas (see for example^^’^^). 

The most classical work in multifluid plasma modeling is due to Braginskii^^, who derived 
fluid equations for a fully ionized plasma, using a Chapman-Enskog closure. Braginskii’s work 
has been successively refined by several authors, with particular emphases on improving the 
transport coefficients and/or including interaction with neutral species^^^^^. Burgers, on the 
other hand, presents a rather general framework for the modeling of elastic collisions^®. These 
include both neutral collisions and charged particle collisions; the methodology is applicable 
for a general system of moment equations beyond the standard five-moment model. Burgers 
also introduces a simplified model for reactive collisions using a Bhatnagar-Gross-Krook 
(BGK) collision operator. 

In the current work, we present a self-consistent model for inelastic collisions within the 
multifluid framework. The model is derived from kinetic theory and obeys the principle of 
detailed balance (DB), which we show to be an essential property to ensure that the system 
approaches the correct equilibrium limit. We focus on characterizing the exchange source 
terms due to collision, namely mass, momentum and energy exchanges (the hydrodynamics 
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and transport flnxes can be added following^^ when considering a non-nniform plasma). We 
will show that in most cases none of these terms can be neglected, and they have complex 
dependencies to microscopic qnantities of the interaction, e.g., mnltiply-differentiated cross 
sections. We present the general description of the collision kinematics and derive the ex¬ 
change terms for the case of excitation and deexcitation processes. Althongh we are mostly 
interested here in electron-impact collisions and atomic transitions, we keep the formnlation 
as general as possible, snch that the application to other species and chemistry (e.g. proton- 
impact, molecnlar vibrational transitions, charge-exchange, etc.) is a straight-forward exten¬ 
sion. The case of ionization and recombination, and other three-body processes, is cnrrently 
nnder examination, nsing the same basic formnlation presented here. 

The rest of the paper is organized as follows. The derivation of the exchange sonrce terms 
is given in Sec. II, by hrst introdncing the description of the transfer integral, and then 
presenting the derivation of the exchange rates in the following snbsections. In Sec. Ill, 
we show the nnmerical evalnation of the mnltiflnid rates, verify the resnits with Monte 
Carlo calcnlations, and perform zero-dimensional calcnlations ntilizing the mnltiflnid rates. 
Finally, conclnsions and a snmmary of the present hndings are given in Sec. IV. We also 
provide several appendices to elaborate on the derivation of the exchange sonrce terms and 
the description of the nnmerical simnlation. 

II. Rate Derivation 
A. Transfer integral 

Let ns consider an inelastic collision between two particles s and t, snch that the particle t 
changes its internal state. The particles s and t are respectively the scattered and target in 
the laboratory frame of reference (LAB). The former will be identihed as the electron and 
the target as the atom, bnt we will keep the general s, t notation nntil explicit assnmptions 
and approximations are made, snch as neglecting terms of the order of the mass ratio ms/rrit 
for hnal expressions. Following Appendix A, the initial velocities are v^, v*, where v = u-|-c 
and u is the flnid mean velocity in the LAB frame, and post-collision valnes are indicated 
by a prime, i.e.: 


s{Vs) + t{vt) s\v's) + t\Vt) 


( 1 ) 


4 


We make here two assumptions: 1) the collision produces only two particles, which may or 
may not belong to the same fluids as the initial reactants, and 2) the masses of individual 
particles are the same before and after the collision, e.g. m'^ = rUg, such that mass conserva¬ 
tion is automatically obtained. Both of these will be revisited in a follow-on paper, dealing 
with ionization and recombination. Defining the energy transfer to and from the internal 
modes to be represented by Ae, we have the following energy conservation constraint on the 
relative velocity g where g = (see Appendix A): 


g' = g'" + 


2A£ 


( 2 ) 


For excitation, the transferred energy is a positive and fixed value Ae = £*, the energy 
gap between the levels, while for ionization it is a continuum of values: Ae G [£*,£], where 
e = is the available kinetic energy in the center-of-mass (COM) frame. In the limit 
Ae —)■ 0, the collision is elastic. We will keep the same relations for the reverse process, for 
which the primed variables are post-collision and non-primed refer to pre-collision, such that 
for deexcitation, Ae = —e*. 


We can then define a transfer integral of the collision operator between the two species s 
and 


= HgTit / (fvgd^vt fsft g duj{vg, vt; v'^, v'*) 


( 3 ) 


where g is the magnitude of the relative velocity {g = |g|), du is the differential cross sec¬ 
tion (DCS), and '0 is any moment variable exchanged during the collision. We now follow 
Appendix B, starting with the following transformations: 


V* 

g 

and 7 


V - U + 7g 

g - w 

KTt-Tg) 

mgTt+mtTg 


,, ^ MTgTt ^2 ^ 2kT* 

mgTt+mtTg ^ M 

~ mgTt+mtTg ^ 2kf 

T =- TV - a =- 

M a 


(4a) 

(4b) 

(4c) 


where the relative mean velocity w = — u*. The product of the two Maxwellian distribu¬ 

tions fg ■ ft is expressed in terms of the product of two other Maxwellians, fv* ■ fg, for the 
COM velocity and relative velocity respectively. Inserting (B.14-B.16) into (3), the transfer 
integral can be written as follows: 


= ngUf 


— S 

9 nA 


d^V* 




TTza 


9 n/A 


d^ge / i^duig] g') 


TT2a 


( 5 ) 
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Figure 1. Frame rotation and relative orientation of (a) w and g and (b) g and g'. The rotation 
operator matrix R{'~p, 9) (or i?(p, x)) is defined such that g = R{(p, 9) • w and g' = i?(p, x) ' g- 


Note that in the COM reference frame, the DCS only depends on the relative velocities, i.e., 
dui^Vg, v^; Vg, v'i) = du{g] g'), and can be expressed as: 

du{g]g') = agt{g,Q')dQ' ( 6 ) 


where D' is the solid angle between the initial and final relative velocities, i.e., dVL' = dpdcosx 
with g-g' = gg' cos x- Without loss of generality, we can now choose a reference frame (LAB) 
such that the relative mean velocity w is aligned with the z axis, as shown in Figure 1. Thus, 
the unit vectors g, g' are obtained by subsequent rotations of the (x, y, z) frame. Using the 
abbreviated notation = cos(p, = sin(p, etc, we define this rotation operator by the 
matrix: 




R{^,9) = 


-se 0 ce 


and g = — = 6) ■ z = 

9 


^ C^Sg ^ 

^ip^O 

\ c, J 


( 7 ) 


Similarly, the post-collision relative velocity is rotated by the angles (p, y), such that g'= 

R{p,x) ■ g- 


Using d^g = g'^dgdipdcg, and equation (6), the transfer integral can be written as: 

■ j d^Y*fv* ■ j dgg^e-^R^"- 

dpdce f dpdcy^pagt{g,id') 


712 a'^ 


( 8 ) 
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Let us now assume that the moment variable can be expanded in terms of powers of V*; 


^ = a + hY* + cV*2 + ... 


where a,b,c... are functions of the remaining velocity variables, and let us perform the 
integration over V*. Note that we have: 


d^V*fv* = 1 


£V*V*fy. = 0 


the latter by reasons of symmetry. Thus, as long as V’ does not contain terms quadratic (or 
higher) in V*, a condition satisfied throughout this work, we can eliminate the integration 
over V*, keeping only the terms which are independent of V*. Also by symmetry, the DCS 
ast does not depend on the angle p, and we can write: 


ct«( 9>S2') = ■ 6(j,\') s.t, dp dc^g(g,x) = I 


( 9 ) 


where Q is the angular-dependent DCS. More generally, we will define the averaging of any 
function -0 over the scattering angles as: 


"+i 


W')n, =2!r/ *,V’6(9,\') 


( 10 ) 




A trivial integration over ip yields: 






( 11 ) 


We now define the following, normalized energy variables. 


^ = 




A = 


1 2 


z = max 0, 


Ae 


( 12 ) 


kT kT \ kT, 

where T is defined in Appendix B. Using g^dg = 2ede/p^ and a further change of variables, 
we finally obtain: 


, 8kT 

= nsTit I - e 



0 + 1 


dzze-‘W„{z) - -j . {^)„, 


(13) 


with gf a thermal velocity based on the average temperature T. Note that we have left the 
variable i/i undetermined, and since it could potentially depend on all integration variables 
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{z,x,d), must be kept inside all integrals. We will see what simplifications can be made 
next, depending on which moment variables we are integrating. 

The lower limit of integration, z*, is zero for elastic collisions or for exothermic reactions. 
Thus, equation (13) is a general formula, which applies equally to excitation {z* > 0) and 
deexcitation [z* = 0)^®. Let us first consider an excitation collision; 

s(vs) + t{Ei, Vi) s(v'J + t{Eu, v'i) 

where £ and u denote the lower and upper energy states, respectively. We made here the 
assumption that both states (£, u) belong to the same fluid, so that the particle indices (s, t) 
are kept the same, but this is not necessary. From eq. (2), energy conservation implies 
Ae = Eu—Ei > 0. We can then define normalized energy variables for this case, x, x' and x*\ 

X* = >0; x' = X — X* > ^ (14) 

kT 

Note that x {x') is the normalized kinetic energy of the initial (final) products of excitation 
respectively, and that x* is the normalized energy threshold, always a positive quantity. For 
excitation, we can use (13) with the following identifications; 

z ^ X, z* = X* , z' = x', and rit = rii (15) 

For a deexcitation collision, we have the reverse (u ^ £), i.e. Ae < 0; 

s(v5) + t{Eu, vt) -)■ s(v^) + t{Ei, v't) ( 16 ) 

Equation (13) is again still valid, if we now make the following identifications; 

z = x', z* = 0, z' = X, and Ut = Uu (17) 

Therefore, in all cases the variable x always refers to the larger kinetic energy (before ex¬ 
citation or after deexcitation) and x' refers to the smaller value (after excitation or before 
deexcitation). We can therefore define two cases of (13); 

poo 1 ^ + 1 

'^l^ = nsnigfe~^ dxxe~'"wl^{x)-- • (-0)^, (18a) 

Jx* 2 

= ^snugfe~^ dx'x'e~^'ai^{x') ■ ^ j (18b) 


where the superscripts i indicate excitation and deexcitation respectively (note the change 
of subscript from st to si for excitation, and su for deexcitation). It is worth pointing out that 
the averaging over the scattering angle, i.e, has to be done with the corresponding 

angular-dependent DCS Q, e.g., = 27t J for excitation. However, from time 

reversal we have Gl^ = Gju^ simplicity, we do not differentiate {'ijj) between excitation 

and deexcitation. Note that the integration over the 6 angle remains to be performed. If 
the moment variable '0 can be expanded in terms of power of cos6', we can then define the 
following set of functions: 

C^"^(0=-A4 rdyy^e^^y (19) 


'-1 


where A4 is a normalizing factor. In particular, we have: 
1 


C'“>K) = 


/-I 25 

'• + 1 o 

.2Cy _ 




cosh (2,^) — 


sinh(2,^) 

2e 


s.t.: hmC*'°^ = l 

(20a) 


s.t.: lim = 1 

(20b) 


In the CR model, each internal state is treated as a pseudo-species, so the rate of change in 
number density for each state {ng, riu) is taken into account separately. We can now examine 
the specific form taken by the transfer integral, according to the chosen moment variable, 
starting from (8), (13), or (18). 


B. Zero*^-order moment: number density 

The rate of change of the number density due to an inelastic collision of type (1) can be 
obtained by setting -0 = 1 in (13), so the average over all the scattering angle is trivially 
removed: 

/ oo P+l 

dzze~^ast{z) ■- J dcee^^""^ ( 21 ) 

The integration over dcg yields the function defined in eq. (20). We can now express the 
rates for transitions between two atomic levels £, u, by making the appropriate substitutions 
for the energy variables. For the case of an excitation (£—)■«), and according to eq. (15), 
we define the variable x as the normalized kinetic energy of the reactants (s, t) in the COM 
frame, prior to the collision: therefore in this case, z = x, z* = x* > 0 and Ut = ng. Thus, 

/•OO 

= nsnigfe~^ / dx x e~^ al^{x)C‘^^\\^Xx) (22) 

Jx* 
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For deexcitation (m ^ the rate of change of number density follows from (17): 


= UsUngfe ^ / dx'x'e 


Both of these quantities are positive, hence the resultant rates equations are: 


( 23 ) 


dTT'i _ pf _ 


driu 

dt 


and 


dni 

dt 


= +T^ = 


duu 

dt 


In the case of electron-impact processes (s = e), we can neglect terms of order mg/M, and 
for an atomic transition between levels i ^ u, we obtain: 


= HeUiVe e ^ / dx X 6 ^ al^{x) 


(24) 


where Ve = \J In the limit of thermal plasma where multifluid effects are weak, i.e. 
A —)■ 0, we obtain: 

poo 

= rien^Ue / dxxe~^al^{x) (25) 

Jx* 

which is exactly the expected result for a single-fluid plasma. 

Using the Klein-Rosseland relation for detailed balance^®, 

^li.{x)xgt = ai^{x')x'gu (26) 

where are the degeneracies of the lower and upper atomic levels respectively, we can 
write the excitation rate as follows: 


= UsU^gf e ^—e ** 


dx'x'e X{x* +x')) (xiuix') 


One can then easily extract reaction rates, for example: 


Ke = ^si ■ 


pi = -rrjl . 

0/1 I 0/1» 




It is instructive to consider the ratio of these rates: 

-I r/7 1 J^dx' x' e~^'(^^'>{^yX{x'+x*))al^{x') 


'^SU 




(27) 


(28a) 

(28b) 


(29) 


j^dx' x' (\/AU)(jI„(a;') 

The first term in brackets is the traditional Boltzmann equilibrium relation; the second 
term contains the correction due to the multifluid effects, and appears only through the 
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function (20-a). A Taylor expansion near A = 0 yields (with an obvious definition of the 
Boltzmann function B)\ 


VJ 


se 


'^su 


^ -X* 

Ue 

Beu{f) 


J^dx'x'e 


1 + 


2A(a:*+r') 




J^dx'x'e-^' [l + a\^{x') 

, 2Aa;* ,,2^ 

1 + ^— + o(A2) 


(30) 


Thus, we recover the expression for Boltzmann equilibrium in the thermal (single-fluid) limit 
(A —)■ 0). Note that the correction term increases with the energy threshold, i.e. transitions 
between high levels {x* 0) will not be affected very much by the multifluid effects, while 

the impact will be stronger for excitation from low energy levels, with high energy gaps. For 
elastic collisions (a:* = 0), the ratio of rates is exactly given by the ratio of degeneracies. 


C. First-order moment: momentum density 

Consider now the forward reaction (1) and the corresponding loss of momentum to the 
particles with velocity v^. The transfer variable in this case is '4) = ms'Vs^ and starting from 
equation (8), the contribution to the momentum equation is: 

R-7 = ■ [ dgg^e~^^^'^^ast{g) ■ [ (31) 

J J J 

Similarly, the gain in momentum is given by the production of new particles with velocity 

v'.: 


= [ d^Y*fv*- [ dgg^e d"Ia st{g) ■ [ (32) 

J J J 

Using the relation; 

ms(v^ - v'J =/i(g - g') (33) 

we verify that the integrand does not depend on V* and its integration is trivially removed. 
The net rate of change to the momentum density of species s is therefore: 

= [ dgg^e~^^/^^ast{g)- [ dy^dce (g-g%, (34) 

TTZQ;-^ d J 
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Let us consider the last integral over the scattering angle. From Figure 1, the vectors g, g' 
in the rotated frame i], are: 





(cs \ 

s = g g = g- 

0 

; s' = g' ■ g' = g' 



Iv 


V / 


Therefore the integral yields: 

j dn'(g-g')g(g,n')^27rgj (k^g(g,x)g - 2Tig' j dc^c^g(g,x)g 

= [s-s'(cosx)„,] g 


(35) 


(36) 


We must now express the unit vector g in the initial (x, y, z) frame, which is given by (7); 
integration over the (p variable leaves only one component, c^w, yielding: 


= f dgg^e ■ f [g - g'{cos x) (37) 

TTZQc^ J J 

Using (20) to replace the last integral and using the normalized variables (12) leads to: 

Ks =--iJ,wnsntgfe~^ J dz z 2 C^^\\/Xz) (cos x}xj,J (38) 

A similar (but of opposite sign) expression can be obtained for the species of type t, as a 
result of the identity rnaiyt — ^'t) = U-is' — s)- 

We can now specify the type of collision. For an £ u excitation, we follow (15), to yield: 

2 poo ^ 

KI = —-g.w nsn£gfe~^ J dx y'\/x — Vad{cosx)g^,j (39) 


For deexcitation, we follow (17): 


R-t = --jJ^wnsriugfe ^ 


dx' {x') 2 e a\^{x')C,^^\\r^) (\f^-^/x{cosx)xii 


(40) 


Note that the expressions (39-40) are obtained in a frame where w is aligned with the z 
direction and corresponds to the change in momentum density along that direction. Thus, it 
is the component of a force parallel to w, while all components in the transverse directions 
are zero, by reason of symmetry^^. The components in an arbitrary rest-frame must be 
obtained by projecting w. 
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Since the force density is approximately proportional to w, we can group all the other terms 
into the definition of a coefficient, such that 


Rl = -Kl^ius - ut) rJ = +Kl^ius - Ui) 
Ri = - ut) r| = - Ui) 


where and Kj^ are known as resistance coefficients: 




^se = -^t^nsHegfe 
I 2 

Ktu = -^l^nsriugfe 


-A 


dxx^e ^ al^{x) y/x — \fx'{cosx) ^ 

) _ 

dx' {x')^ e~'^' al,i^{x) C^^\V Ax') \fx' — \/x{cosx)^, 


(41a) 

(41b) 

(42a) 

(42b) 


It must be pointed out that when the collision is elastic, i.e., x = x', we recover the expression 
of the momentum transfer cross section often used in transport calculation, i.e., a^(x) = 
a{l — (cosy)) (see for example^^). In the limit of weak divergence of mean fluid velocities 
(A —)■ 0) and isotropic scattering {Q{x) = l/4vr), we have: 


/ dxx'^ali{x)e' 


(43) 


Again, using the Klein-Rosseland relation, the excitation resistance coefficient can be written 
as: 




~ 1 2 

^iu{T) -iinsnigfe~^ / 

-I o Jo 


dx'x'x’^e y/x — {cosx) 


(44) 


As in the case of the zero-th order moment, we define the momentum exchange rates by: 

kI^ = pinsUiKlf^ (45a) 


Ki. = 


finsHuK; 


(45b) 


In the case of weak divergence of mean fluid velocities and isotropic scattering, the ratio of 
the rate coefficients for the forward and backward processes is approximately: 


K 


si 
1 

f^SU 


Beuif) 


dx'e ^'x'{x* + x') [l + |A(x* + a;')] (tI^{x') 


(46) 


dx'e~^'x'‘^ [l + |Ax'] al^{x') 

Note that there is an additional contribution from high-order moment from the expansion. 
This can be seen by further expanding the integrand of the numerator: 


x'{x*+x') 


1 -h -A(a;*+x') 
5 


J2 


= X 


1 + -Xx' 


I * f 

+ x x 


1 + -A(a:* + 2x') 
5 
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such that 




4 ^'x' [1 + |A(a;* + 2x')] (t\^{x') 



dx'e~^'x'‘^ [l + fAa;'] d'l^{x') 


( 47 ) 


Note that even as A —>■ 0, the correction term does not vanish. Therefore, there is no 
equivalence between the resistance coefficients of the forward and backward processes in 
the limit A ^ 0. However, this is perfectly understandable; note that the correction is 
proportional to the energy threshold, and since kinetic energy must be removed from particle 
s in order to achieve excitation, but not for deexcitation, there must also be an imbalance 
in the momentum exchange rate. As expected, this imbalance vanishes for elastic collisions 
[x* 0), and the rates are consistent with the detailed balance of the mass exchange. In 

all cases, detailed balance is enforced through relation (26) at the microscopic level. 


D. Second-order moment: total energy density 

The net rate of change of total energy of species s can be obtained by setting 'ijj = 
^nis (v'^ —into equation (11): 

■ J J dg ast{g) ■ J (48) 

Using the transformation defined in Appendix B, 

= h (g'-g) • [V* + U- 7 (g-w)] - ^Ae (49) 

The integration of the hrst term in the square bracket is zero since J d^V* V* fv* = 0. One 
can easily see that the second term in brackets is simply ■ U by comparing with (34). 
Similarly, the last term in (49) is identihed as —{mt/M)Aerst- The third term in brackets 
involves the following dot product: 

(g'-g) ■ (g-w) = g ■ g' - w ■ g' - 5(2 + w ■ g 

= gg' cos X — wg'{g' ■ w) — g'^ + wg cos 9 (50) 


We can now perform the averaging over the scattering angle; in particular, we have 


jdVt'Q{x) g'-w= (cosy)^, cosd 


14 








so that 


j (g'-g)-(g-w) = {g'{cosx)n,- 9 ){ 9 -wcose) (51) 

Thus, we obtain, after integration over ip: 

Q, = U-R,-^A£r,i-/i7^e-^ [dg Te (52) 

ivi J 

where the last angular integral is 

Ie = ^ jdc0e^^^‘'^^^\g'{c^)^,-9){g-wc0) 

= (g'ic^).,-9) {4 /- 11 dc0C0e^^-^^/^^ | (53) 

= i9'{c^)n'-9)9 


Therefore, after the change of variables g ^ z: 
Qs =U • Rs — — AeTsi 


+ fi 7 ^^|^e jdzzie "a„(z] (yz-\/7{cos (^C*°’(^V)-yC*‘’(\/V)^ 


(54) 


The factor in front of the integral can be re-arranged to yield; 


ynsUtgfe ^{2kT) 


Using the identity y{2kT) = ■^k{Tt — T*), the hnal result has a traditional form; 

Q, = U • R, - ^AeT,* + - Ts) 

with the thermal resistance coefficient dehned as; 

Ja = n,ntgte~’' jdzzie~"a,i{z) (^Vz-\/z'{cosx)„,^ ('c*”*(VA2) - 


(55) 


(56) 


This result is general, and we can now make the usual substitutions for excitation; 


Q: = U . R; - ^£'7, + Jl,‘^k(T, - T, 


(57a) 


= UsUegfe ^ jdxx^e (^x/x-V^{cosx)xi,^ j (57b) 
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In the case of deexcitation, we can still use the general formula (54), except that in this case, 
Ae = —e*. This can be seen if we start from eq. (48), which gives us: 


J 

+ /i 7 ^e-" [ • (g-w) (58) 


vrao;'’ 

^ M 


Again, one can easily recognize the standard formulae: 

Qi = U ■ Rt + + JL^m - Ts) (59a) 

Jju = UsHugfe-^ Jdx'{x')h-^'ai^{x') (y^-^/x{cosx)^,^ 

(59b) 

which could also be obtained directly from (55 - 56), with the usual substitutions (17). 


We can also express the source term for particle t, using: 


1 / /2 2\ ^ f /-tT /\2 /’i;- \2 


rus 

2M 


Ks' -g ) -/xV- (g'-g) 


(60) 


m. 


= -/^(g'-g) ■ [V* + U - 7 (g-w)] - —Ae 


Comparing with (49), we easily obtain (note the inversion of Tg and Tt in the last term): 

o; = U . R; - + Jl, ^k{T. - T.) (61a) 

Qt = U . Ri + + Ji ^k(T. - T.) (61b) 

with ^ = —RA h Combining both s and t fluids, the only term remaining is the loss of 
energy equal to the energy gap between the levels, as expected. Note also that this energy 
loss is distributed to the respective fluids according to the ratio of masses, such that the 
lighter element receives the major contribution. This is also an expected result, similar to 
the energy exchange due to elastic collisions, and due to the kinematics of collision. 


In the limits of near single-fluid (A —)■ 0) and isotropic scattering, we have: 


2AA 


~ UsUigf ( 1 - yj y dxx^ crL(a^)e 


(62) 
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The thermal relaxation rates can be extracted similarly: 


(63a) 

(63b) 


JIi = risna]^ 

JL = nsUujt 


The ratio of the thermal relaxation rates can be written as: 

^ dx' e~^' x' [l — |A + |A(x* + 2x')\ 


■t 

jsu 




Beuif) 



l + x* 


(7^ 
^ su 


(64) 


dx' e-‘' 1 “ [1 - |A + |Ax'] ffJ„ 

Similar to the case of momentum transfer rates, the correction terms do not vanish when 
A —)■ 0 due to contribution from high-order moments. 


III. Numerical Results 

In the following sections, we carry out a numerical evaluation and verification of the exchange 
rates derived in IIB, IIC and IID, for the case of free electrons interacting with hydrogen 
atoms; these processes include electron-neutral elastic collision and electron-impact exci¬ 
tation and deexcitation. Ionization and recombination are currently omitted and will be 
included in future work. For comparison purpose, we also show the exchange rates due 
to Coulomb collision, i.e., electron-hydrogen ion, which is the dominating elastic exchange 
mechanism for plasma with high ionization fraction. 

The notations are slightly modified to better distinguish each type of interaction. We use 
superscripts (en), (ei) and {xd) to denote electron-neutral, electron-ion (Coulomb), and 
excitation/deexcitation (as a whole) collisions, respectively. These grouped notations are 
useful, for example, when looking at the net momentum (or energy) transfer due to each 
type of collision. The symbols I are still retained to indicate individual excitation and 
deexcitation rates. For each transition between two atomic states, we use the convention 
of indexing the final state on the left, and the initial state on the right, i.e., {f\i)- For 
example, is the forward excitation rate from i to u, and is the reverse process. 

The energy levels and cross sections models for atomic hydrogen are given in classical form 
and summarized in Appendix C. 


A. Reaction Rates 

All the exchange rates (mass, momentum, energy) can be tabulated as a function of two pa¬ 
rameters: the average thermal temperature T, defined in appendix B, and a non-dimensional 
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parameter A, defined in eq. (12), which corresponds to the relative mean kinetic energy. For 
convenience, we also dehne an equivalent drift temperature = AT^^, such that all the rates 
can be tabulated in terms of two temperatures. Since the mass ratio between the electron 
and the atom is very small {rUe <C Mh), we can drop terms of order rUe/M and arrive at 
the following approximations: fx ~ rrie, T ~ Tg, ~ ATg and ~ Ug = There¬ 

fore, all the exchange rates for electron-induced collisions (both elastic and inelastic) can be 
numerically evaluated in terms of the electron temperature Tg and the drift temperature T^. 

Figure 2 shows example calculations of the zeroth-order reaction rates, dehned in eqs. (22), 
(23) and (28), for electron-impact excitation and deexcitation between the hrst three atomic 
states of hydrogen. These rates exhibit a similar trend for the range of temperatures plotted 
here, that is, starting from low temperature, the rates hrst increase, reaching a plateau and 
then decrease as temperature further increases. The value at which the rate is maximum 
is very close to the threshold temperature of the transition. This trend holds both in the 
direction of increasing thermal Tg (x-axis) or drift temperatures Ty, (y-axis). 

It is clearly shown from Figure 2 that the reaction rates can be signihcantly different from the 
thermal limit when the relative mean velocity between two huids is signihcant. In particular, 
one sees an increase of the reaction rate in the low temperature regime where Tg and Ty, are 
small compared to the excitation temperature of the collision; this enhancement corresponds 
directly to the form of the cross sections. Therefore, one can expect that signihcant deviation 
from the thermal rate occurs when the mean kinetic energy is of the same order as the 
excitation temperature. This indicates that excitation and deexcitation among low energy 
states with large threshold energies are more sensitive to the multihuid effects. This is 
consistent with the prior statement we made when examining the ratio of forward and 
backward rates. 

Figure 3 shows the forward and backward reaction rates of the hrst transition between the 
ground state and the hrst excited state as a function of the thermal temperature for several 
drift temperature values. Even at very low thermal temperature (Tg ~ 0.1 eV), one can have 
signihcant excitation {w ~ 10“^“^ m^/s) due to a high drift temperature. In addition. Figure 3 
also shows that in the limit A —?■ 0, the multihuid rate, as formulated here, converges to the 
expression for thermal limit, given in eq. (C.3), as expected. Figure 4 shows the reaction 
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Figure 2. Multifluid reaction rates for electron-impact excitation/deexcitation collisions as a func¬ 
tion of two temperatures. For better color contrast, all the values smaller than minimum value of 
the colormap are white out. 
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(a) Excitation (1 —2) (b) Deexcitation (1 ■(— 2) 


Figure 3. Multifluid reaction rates for electron-impact excitation/deexcitation collisions: lines with 
symbols correspond to different values of the drift temperature. The solid line corresponds to the 
exact solution for the thermal limit given by equation (C.3). 

rates as a function of the drift temperature for several thermal temperature values; here 
we can identify a different asymptotic limit of the rate. In the limit A —)■ oo, the electron 
velocity distribution function approaches the form of a delta function centered at the relative 
mean velocity w, and the reaction rates approach the beam limit given by eq. (C.5). It must 
be noted that in the numerical integration of the multifluid rates, e.g., eq. (22), the energy 
grid X needs to be rehned near the value of the the mean kinetic energy A to avoid numerical 
error due to the integration over a delta function. 

B. Momentum and Energy Exchange Rates 

We now compute the momentum and energy exchange rates due to both excitation and 
deexcitation, and compare with those due to elastic collisions (electron-neutral and electron- 
ion). Recall from eq. (55) that the total energy transfer include three terms: the hrst term 
due to work done by friction in the COM reference frame, the second term due to thermal 
resistance (or thermal relaxation), and the last one due to heat release/absorption due to 
chemical reaction. There are three different rate coefficients associated with each of these 
processes, namely the momentum exchange rate k, thermal relaxation rate j, and the reaction 
rate vj. In this section, we will focus on examining the momentum exchange and thermal 
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(a) Excitation (1 —2) (b) Deexcitation (1 2) 


Figure 4. Multifluid reaction rates for electron-impact excitation/deexcitation collisions: lines with 
symbols correspond to different values of the thermal temperature. The solid line corresponds to 
the exact solution for the beam limit given by equation (C.5). 

relaxation rates. The expressions derived in IIC and IID can be readily used for the case 
of elastic collisions by simply setting z* = 0 and z' = z = x. For example, in the case of 
electron-neutral (en) collision, we have: 

(65a) 
(65b) 


Z / 

/ dxx^e-^a^^^\x) [l-(cosx)„,] 


jF”') = tigC ^ / dxx^e [l —(cosx)f,,] 




It can be seen from the previous two equations that in order to compute the rate for the case 
of elastic collisions, we only need the so-called momentum transfer cross section (T^®"'^™'(a;) = 
(1 — (cos x)). These cross sections are available for a wide range of neutral species due to 
their extensive use in calculation of transport properties (see for example^^). The electron- 
neutral collision cross section utilized in this work is taken from Bray and Stelbovics^^. For 
inelastic collisions, we need the full DCS, i.e., both a{x) and Q{x,x) eq. (9) for each 
process. These cross sections are generally not available and analytical approximation is 
needed. For Coulomb collision, we use the analytical DCS from Rutherford’s scattering 
formula with suitable cut-off based on the Debye length^, yielding a momentum transfer 
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Figure 5. Elastic momentum transfer cross section for electrons in atomic hydrogen computed with 
the DCS from equation (68). The symbols are theoretical values from Bray and Stelbovics^^. 


cross section of the form: 


cr 


(ei)m _ 


IbTreo^^ 


In A 


( 66 ) 


where In A is the well-known Coulomb logarithm, 


/y3 \ 1/2 

A = 1.24 X 10^ f 

Typically, In A 5 — 20; for convenience, we take a constant value of InA = 5 for all the 
plots shown here. In this case, the momentum exchange and thermal relaxation rates can 
be obtained in exact forms^*: 


K>“i(r„A) = 


InA 


167reQ(/cTe)2 A 
InA 


f(^) 

Va 


,-A 


TT 


167reg(/cTe)2 
where erf is the typical error function. 


,-A 


TT 


(67a) 

(67b) 


Due to the lack of data of the DCS for inelastic processes, we have used an analytical Born 
scattering approximation for a Coulomb screened potentiaT^, given by the following form; 

C , 1 


Gi£,x) = 


(1 — hcosxy 


h = 


1 + 


21.8 


( 68 ) 
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where C is a normalization constant such that 2ti Q(e, x)dc^ = 1- This angular-dependent 
DCS has been used to compute momentum transfer and thermal relaxation rates for both 
electron-neutral and excitation/deexcitation collisions. We note here that the angular- 
dependent DOS’s Q for excitation and deexcitation have to satisfy detailed balance^*^: 

(69) 

The above condition implies that the probability for a deexcitation collision with an incident 
energy e to have a scattering angle of y is the same as that for an excitation collision with an 
incident energy e+e*] therefore, one cannot independently specify both DOS’s for the forward 
and backward processes. Figure 5 shows a comparison of the computed elastic momentum 
transfer cross sections to the result from a direct close-coupling calculation of Bray and 
Stelbovics^^; the agreement between the two is excellent. Using the angular-dependent DOS 
defined in eq. (68), the exchange rates, e.g., n and j, can be computed for each bound-bound 
transition and summed over all transitions to yield the total rates: 

EE( (70a) 

i u>£ 

^ + 3\f,\u)'^un^ (70b) 

t U>1 

Based on the total frictional and thermal resistance coefficients, we can extract average 
momentum transfer and thermal relaxation rates for all excitation/deexcitation processes as 
follows: 


^ (71a) 

jM) = (71b) 

where is fhe total atomic number density (summation over levels). Note that 

according to our definitions, the exchange rates and contain terms designating the 
population of the excited states, e.g., nk/rin- For comparison purpose, these average rates are 
calculated by assuming a Boltzmann distribution of the atomic states, i.e., Uk = —— 

and Zn = population of the excited 

states is effectively replaced by a Boltzmann distribution characterized by a temperature 
Tb- This step is only done for the comparison shown below. In a CR calculation, the 
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detailed population of the atomic states (equilibrium or not) is known and the exchange 
rates are computed for each transition as specihed in eq. (70). 

Figure 6 shows a comparison between the momentum exchange and thermal relaxation 
rates for three different processes: electron-neutral, electron-ion (Coulomb), and excita¬ 
tion/deexcitation collisions. For clarity, the rates due to Coulomb collision are only shown 
for value of Tg > 1 eV. Two important observations can be deduced from this plot. Firstly, 
when the atoms are cold, i.e., Tg is low, the inelastic exchange rates are much smaller than 
elastic ones. This is due to the fact that when Tg is low, only collisions between the ground 
state and a hrst few excited states are signihcant; the rates for these transitions are low com¬ 
pared to the others due to larger energy threshold. As the atoms are being excited and Tg 
increases, transitions among highly excited states become signihcant, leading to an overall 
increase in the total rates. These rates eventually exceed those due to elastic collisions as can 
be noticed in region (I) in Figure 6 for the dash-dotted and dotted lines. Secondly, the rates 
due to inelastic processes tends to have a slower drop-off at high temperature compared to 
the elastic rates, which suggests that at sufficiently high temperature, the main momentum 
and energy transfer mechanisms (region (II) in Figure 6) are due to inelastic collisions. In 
addition, one can also observe that the thermal relaxation rates decrease at low thermal and 
high drift temperatures (red solid curve with and without circles). This is due the cancel¬ 
lation of two terms inside the last bracket of eq. (56); the second term has a multiplicative 
factor of A hence its magnitude is increased at high drift temperature. Thus, we are led 
to the important conclusion that one cannot neglect momentum and energy transfer due to 
inelastic collisions. The only justihcation for neglecting these terms is when the atoms are 
cold and thermal temperature is low, but both of these conditions will not be realized in 
most practical systems. 

C. Verification 

The accuracy of the derived formulas of the exchange source terms are verihed against direct 
evaluation of the full transfer integral (3) over six dimensional space using Monte Carlo 
method. The procedure for the Monte Carlo integration is as follows: (1) sample different 
pair of particles (one atom and one electron) from two different Maxwellian distributions, 
(2) compute the exchange rate due to each sample pair, and (3) accumulate these rates. The 
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T, [eV] 


TAeV] 


(a) Momentum exchange 


(b) Thermal relaxation 


Figure 6. Comparison of the momentum exchange rates n and thermal relaxation rates j due to 
elastic and inelastic collisions (excitation and deexcitation). Different colors of the line indicate 
different values of the drift temperatures in eV: = 0.1 (blue), 1 (green) and 10 (red) eV. Solid 

lines indicate rates due to e“-H elastic collisions, solid lines with circles indicate e“-H'’' (Coulomb) 
collisions, and three sets of broken lines are used for the inelastic collisions. The rates for inelastic 
collisions are defined in eq. (71) and are: dashed, dash-dotted and dotted lines for the case of 
Tb = 0.5,0.7 and 0.8 eV, respectively. 

sum of these rates will follow the correct probability distribution function of the samples. 
For brevity, we only show an example calculation of the zeroth-order reaction rate for an 
excitation and deexcitation from levels 2 and 5. Figure 7 shows a detailed comparison 
between the rate expressions obtained from eqs. (22) and (23) against Monte Carlo results 
with an excellent agreement. Similar agreement is obtained with other sampled transitions, 
giving us complete conhdence in the accuracy of our derivation of the multifluid rates from 
kinetic theory. 

D. Zero-dimensional Calculations 

We conducted zero-dimensional (OD) calculations for a constant-volume (isochoric) system 
to support our hndings in sections 111 A and lllB, which include the following; (1) reaction 
rates can be enhanced when the relative mean drift velocity is signihcant, and (2) momentum 
transfer and thermal relaxation due to inelastic collisions are non-negligible. The system 
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(a) Excitation (2 5) (b) Deexcitation (2 t— 5) 

Figure 7. Comparison of zeroth-order reaction rates with Monte Carlo integration of the full transfer 
integral. Lines are from our derived expressions. Symbols are Monte Carlo results. For each line, 
the ratio of the drift and thermal temperature A is fixed and the values are shown in the figure. 

contains hydrogen atoms and a small fraction of free electrons and ions. Here we consider an 
electron-hydrogen two-fluid system. The relaxation between the heavy particles are assumed 
to be inhnitely fast such that they can be described as a single heavy fluid with a bulk velocity 
u/j and a temperature T^. The governing equations are described in appendix D, and the 
resultant system of ordinary differential equations is solved using the RadauS method of 
Hairer and Wanner^®, which is ideally suited for stiff problems. 

The initial conditions for the number densities and temperatures of the heavy particles and 
free electrons are summarized in Table I. Initially, all the atoms are at the ground state 
(denoted by A; = 1 in table I) with a translational temperature of 0.3 eV. The atoms and 
ions are assumed to be at rest, i.e., their mean velocity is zero. The free electrons have a 
temperature of 2 eV, and their mean velocity is varied to demonstrate the multifluid effects. 
While the electron number density is fixed for all the test cases shown here, the ion density 
is varied to examine the role of Coulomb collisions. In all the test cases, we include the hrst 
10 atomic levels of hydrogen according to the model described in Appendix C. 

In the hrst test, we perform the calculation with various initial mean velocities of the elec¬ 
trons, or equivalently, the drift temperatures T^. The number density of ions is set to zero 
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number density 

temperature 

atomic 

Uk = 0.9ni for k = 1 

0.3 eV 

states 

Uk = otherwise 


ion 

ni = fue] f = 0,0.2,1 

0.3 eV 

electron 

Ue = O.lrit 

2 eV 


Table 1. Initial conditions of OD test cases. The ion density is varied for each test case. For all cases, 
the total atomic density rit is 10^^ m“^. 


in this test, so there is no Coulomb collision. Figure 8 shows the time evolution of the 
number densities of the excited states during the isochoric heating process for two extreme 
cases: = 0.01 and 10 eV; the former corresponds to a single-fluid calculation of a thermal 

bath with a warm electron population, and the latter corresponds to a situation where an 
electron beam is injected to the system. One can clearly see from Figure 8 that there is an 
enhancement to the excitation process, indicated by an early increase in the population of 
excited states, due to the presence of a non-zero mean velocity of the electrons. The same 
argument can be made from Figure 9, which shows the time evolution of the temperatures for 
three cases of different initial T^. In this plot, the Boltzmann temperature Tb indicates the 
degree of excitation of the atom. It must be pointed out that the enhancement in excitation, 
however, persists on the time scale of the momentum relaxation process, which is indicated 
by a drop in at approximately 3 x 10“^ sec as shown in the bottom plot of Figure 9. After 
this time, the momentum of the electrons is completely absorbed by the atom, signifying a 
change to single-fluid kinetics. In all the test cases, the excitation proceeds at a time scale 
much smaller than the resolution of the hgures, i.e., Tb approximately goes from 0 to 0.7 
eV in 10“® sec. As mentioned before, when Tb is sufficiently large, the momentum exchange 
and thermal relaxation rates from inelastic collisions cannot be neglected. 

In the second test, we specihcally identify the effect of inelastic collisions. Figure 10 shows 
a comparison of the temperature evolution for three test cases, all with a very small initial 
drift temperature T^ = 0.01; the solid lines correspond the solution with both elastic^® 
and inelastic exchanges, and the dashed lines to the solution without inelastic exchanges. 
In these test cases, we slowly introduce the ions which effectively increase thermalization 


27 











Figure 8. Number density of excited states during a zero-dimensional chemistry test. Solid lines 
correspond to the solutions where the relative mean velocity is very small (T^, = 0.01 eV). Dashed 
lines correspond to the solutions with a large relative mean velocity (T^ = 10 eV). 
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Figure 9. Time evolution of the temperatures for several test cases with different initial drift 
velocities: solid line (T^ = 0.01 eV), dotted line {Ty^ = 1 eV), and dashed line (T^ = 10 eV). The 
bottom plot shows the evolution of the drift temperature Ty,. 
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between the heavy particles and the electron via Coulomb collisions. The friction is negligible 
since ~ 0, and there are only thermal relaxation and heat release from reaction. It is 
clearly shown that inelastic collisions do contribute to the total energy transfer between the 
electrons and atoms, leading to a faster temperature equilibration T/^ —Tg —Tg. These effects 
are most noticeable when the ion density is low due to a weaker contribution of Coulomb 
collisions. Nevertheless, we observe under-prediction of temperature relaxation in all cases 
when inelastic energy exchanges are omitted. It is interesting to point out here that while 
Te — Tb equilibration is due to the heat release/absorption term, Tg — Th equilibration is 
only due to the thermal relaxation term. This result is quite intriguing, since the impact of 
inelastic collisions on thermal relaxation is normally neglected in most of single-fluid multi¬ 
temperature calculations. 

Figure 11 shows the results for the same three test cases above but now with an initial 
drift temperature = 10 eV. One can make the same argument that inelastic collisions 
further enhance the temperature equilibration process. However, it is worthwhile to point 
out that during the time period 10“® < t < 10“^, the electrons, although being decelerated 
due to friction, are also heated due to the work done by the same force (first term on the 
right hand side of eq. (D.5a)). This thermal heating term contains contribution from both 
elastic and inelastic collisions, which explains why the electron temperature drops faster 
towards equilibrium when inelastic exchanges are neglected. Finally we should emphasize 
that the thermal equilibrium between all components (Te,Te,T/i) can be obtained purely 
from inelastic collisions, as a result of properly accounting for detailed balance in our model; 
elastic collisions are not required to achieve equilibrium, but of course are needed to obtain 
the correct rate of relaxation. 

IV. Concluding Remarks 

We have presented a model for inelastic collisions for electronic excitation and deexcitation 
within the context of a multifluid description of a plasma. The model is rigorously derived 
from kinetic theory and is applicable to any multifluid plasma, irrespective of the mass 
ratio, and strictly obeys the detailed balance principle. The appropriate mass transfer rates 
and momentum and thermal resistance coefficients are derived, and are found to satisfy the 
proper asymptotic limits, such that in the limit of vanishing energy gap, the well-known 
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Figure 10. Time evolution of the temperatures for the case with (solid) and without (dashed) the 
exchange terms (momentum and energy) due to inelastic collisions. The ion number densities are 
specified for each case. Initially, = 0.01 eV. 


expressions for elastic collisions are recovered. 

Nnmerical evalnations of the mnltiflnid rates are carried ont for a two-flnid electron-hydrogen 
plasma nsing Bohr model for the energy levels and semi-classical cross sections. Several nn¬ 
merical tests were performed in a virtnal (zero-dimensional) test cell, and both the known 
thermal and beam limits were correctly recovered from the model. We also fonnd that in 
some plasma conditions of interest, the contribntion of the inelastic collisions to the resis¬ 
tance coefficients is signihcant, contrary to the nsnal assnmptions made in cnrrent mnltiflnid 
models. While this work is focnsed on two-body reactions, extension to ionization and re¬ 
combination processes is cnrrently nnder examination, and will be presented in a fntnre 
article. 
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Figure 11. Time evolution of the temperatures for the case with (solid) and without (dashed) the 
exchange terms (momentum and energy) due to inelastic collisions. The ion number densities are 
specified for each case. Initially, = 10 eV. 
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potential is lowered as a result of interaction with the plasma (Debye shielding) and quan¬ 
tum uncertainty. In practice, the truncation is accomplished at a lower limit still; for the 
current purpose, details of this truncation procedure can be ignored. Suffice to say that 
the series extends to a number n = M, which can be considered large, e.g. 0(100). 

Appendix A Collision kinematics 

Let us consider an inelastic collision between two particles s and t, such that the post-collision 
particles can have modified internal states. The process is formally described as the relation 

s(v*) -F t{wt) -t s'(v'J -h t'(v't) (A.l) 

Note that only two particles are produced by the collision. The initial velocities are Vs,vt. 
The mean fluid velocity is u, such that u = (v) = J dfwf{\) and a thermal velocity 
c = V — u. By definition, we also have (c) = 0. 

The collision can be transformed to the center of mass (COM) reference frame, moving with 
velocity V with respect to the LAB frame. Similarly, we can also define a mean velocity 
of this COM frame as U. The subsequent Galilean transformations yield the following 
definitions: 


V = 

ms^s + mtYt 

g = ^s-^t 

(A. 2 a) 

M 

u = 

msUs -f mtVLt 

W = Us - Ut 

(A. 2 b) 

M 
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where M = ms+mt. The inverse transformation yields: 


= V + —g u, = U + —w (A.3a) 

Vi = V - Ui = U - (A.3b) 

Mass conservation imposes the relation ms+rnt = M = m'g + m[. For the case of two-body 
processes such as excitation of internal states, the masses are individually conserved, i.e. 
m'g = ms-iTn[ = mt. Expressed in the COM frame, momentum and energy conservation yield, 
respectively: 


MV =MV' (A.4a) 

^MV^ + ^/xg2 =^MV'^ + ^/ig'^ + Ae (A.4b) 

where fi = msmt/M. Therefore, we have the following constraints: 

9 Ar 

V = V' and g^ = g'^ -I - (A.5) 

/i 

For an excitation between two atomic levels, the transferred energy is a fixed value Ae = 
e* > 0, the energy gap between levels. For a deexcitation, we use the same formulation as 
above (i.e. post-collision variables indicated by a prime), but in this case, Ae = —e* < 0. In 
the limit Ae —)■ 0, the collision is elastic. 


Appendix B Separation of variables 


Consider the Maxwellian velocity distribution functions (VDF) of each particle type, nor¬ 
malized to unity, e.g. (recall that c = v —u): 


/s(v,) 


/ V r 


(B.l) 


and similarly for ft. The averaging over initial states will yield a product of these two 
distributions: 


/s(v.)/i(vt) = 


m.. 


27r/cT, 




(B.2) 


where the argument of the exponential function is, from (A. 3): 


A = 


m.o 


V-U + —(g-w) 
2kT, L M ^ 


+ 


m 

2kTt L 


V-U 


ms 

M 


(g-w) 


(B.3) 
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Following Burgers^®, this expression can be simplified with an appropriate transformation of 
variables; since the basic procedure will be used elsewhere, we describe it below. First, we 
define the following variables 






^ 2kTp 


g = g - w 


(B.4) 


such that 


rrit . 


A = ^ (V-U) + ^gJ +/?,^(V-U)-^g 


m.. 


= (/3,+A)(V-Uf + 


2 2 
mi ^ mz 


r+ 2 


rrit nis 




Let us define: 


V* = V-U + 7g 


(B. 6 ) 


and comparing the expression 


= Ws + my-w + (/9.+A)yr + 27 (/?.+a)(v-u) ■ g (b.?) 


with (B.5), we can choose the appropriate value of the coefficient 7 to eliminate the dot 


product from A: 



We then obtain complete separation of variables: 


(B. 8 ) 


^=(/3,+A)V*' + 




M2 ^‘M2 


1 

/^s +A 




The term in brackets is easily simplified: 


/3s/3t 
/3 s + A 


(B.9) 


(B.IO) 


We can now define effective, average temperatures: 

_ rus _ M nisTt + rritTs _ M 

^ MTsTt ~2kf* 
/3s/3t _ M TsTt ^ /i 

A + A 2fc TsTt msTtTmtTs 2kf 

and 7 becomes: 

_ H Tt — Tg _ Tf — Ts 

^ ~ M f ~ ^ msTt + mtTs 


(B.lla) 

(B.llb) 


(B.12) 
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To summarize, we have performed the following change of variables: 


V* 


V-U + /i 




msTt+rutTs 


M- 


TsTt 


msTt+mtTs 


g = g - w 

^ _ msTt+mtTs 

M 


(B.13a) 

(B.13b) 


These are the same expressions found in^* (pp. 45-46) (with an occasional change of naming 
convention) for which it is easy to verify that the Jacobian of the transformations is unity, 

i.e. 

= d^Vd^g = d^\*d^g (B.14) 


Furthermore, we note that; 


rris 

2kTs 




Wt) 


(/^s+A)^ 




3 

2 



(B.15) 


The product of two distributions can now be written as: 


fs- ft = 


M 


2TTkT* 


exp 




*2 


2kT* 


\27ikTj 


exp 


hg 


2 1 


2kT 


= r(V*)-/(g) (B.16) 


All subsequent expressions can now be simplihed with this separation of variables. For 
example, any operator ^ that depends only on variables expressed in the COM frame (g, g'), 
we have: 

I dVvi/J,^(g,g') = J d^V*nV*)-1 dW(g)^(g,g') (B.17) 

'-V-' 

=1 

Note that this procedure applies equally well for elastic, excitation and deexcitation collisions, 
and that no approximations have been made on the mass ratio. Furthermore, since the 
averaging over initial states only involves the distribution functions for the s and t particles, 
we have not necessarily assumed that the hnal products s' and t' belong to the same fluid 
as the initial particles. 


Appendix C Atomic data and cross section models 

The atomic states of the Hydrogen atom are listed as a function of their principal quantum 
number (n) only, following the Bohr atomic model; the splitting of states with respect to 
orbital and spin numbers is ignored, and all states have a degeneracy = 2n^. The states 
number from n = 1 to cxo and we consider a hnite number of states n = 1 ,..., M < cxo before 
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reaching the ionization liniit^^. In this simplihed model, the energy of each state is given 
as En = Ih (1 — 1/''^^), as measured from the ground state {Ei = 0), and we will denote by 
In = lH (1/n^ —1/M^) the energy required for ionization of level n. 

The classical form of the cross section for energy exchange between a free electron and the 
atom (Hydrogen) is used^. For an excitation collision from level ^ to level u > the cross 
section takes the form: 


= (47ra2)(3ffa) 


{x - xeu) 




(C.l) 


where Oq is the Bohr radius, x is the nondimensional incident energy of the electron, x^u = 
{Eu — Ei)/kTe is the energy gap between I and u, and f^ is the oscillator strength; 

32 1 1 1 


^tu — 


37r\/3 


iw ~ 


(C.2) 


In the thermal (single-fluid) limit (A —?• 0), the reaction rate can be obtained in an exact 
form: 

2 


^lu\e) - (47ra^)(3f&)ne f ^ j 


(C.3) 


where 


\ 2 g 

Vp_ = I -^ , iJeu = - Ei{xiu) and Ei(a;) = 

vrm^ / 




dy 


(C.4) 


X(,u Jx y 

Here, Ve is the mean thermal electron velocity and Ei is the exponential integral. On the 
other hand, in the beam limit (Tg —)■ 0), the reaction rate takes the form; 

2 


= (47ra^)(3h„)ne 


where 


yF (A - Xfa) 

2 XiuX^ 


JJL 

kE 


(C.5) 


(C.6) 


Appendix D Rate equations 

We describe here the governing equations for the zero-dimensional simulation, which describe 
the time evolution of a constant-volume system during a thermochemical relaxation process. 
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Here we made the approximation that rrii nih- The following system of rate 

equations is considered: 


dUe 

dt 

drii 

dt 

drik 


dt 


d 

— (meneUe) 

^{mhrihUh) 

dt 


dEe 

dt 


dEh 

dt 


0 

0 

k 

-(Ue - 

+ (ue - Uh){K^^^^ + 

-u • (Ue - 

‘^k{Te - 
rrih 

+ U ■ (Ue - 

‘^k{Te - Th){J^^^^ + 


where Uh 


rii + J2ken '^k and Ee{h) is the total energy of the fluid: 
3 1 

Ee[h) '^^e{h)kTe{h) T”nie(h)ne(/i)Ue(/i) ' Ue(/i) 


e(h) 


(D.la) 

(D.lb) 

(D.lc) 

(D.ld) 

(D.le) 

(D.lf) 

(D.lg) 

(D.2) 


where £e(h) is defined as the thermal energy of the fluid. The exchange source terms can be 
decomposed into three parts: elastic and inelastic collisions between the electrons and atoms 
and Coulomb collisions between the charged particles. The momentum and energy exchange 
rates for these processes are defined in eqs. (65), (67) and (70). The remaining terms are: 


-pi (^xd'j 

^ k 


E 


-nkUeW 


+E 




(D,3) 


'^{xd) 


u>k £<k 

7^ 7^ ~ 


(D.4) 


i u>£ 


where = E^—Ei. Comparing eqs. (D.ld)-(D.lg) and (D.2), we can also write conservation 
equations for the thermal energies of electrons and atoms: 

dSp 


dt 


= -me(u - Ue) • (Ue - 


2me 

rrih 


k{Te - - 5^ 


xd) 


(D.5a) 
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dsh 

dt 


(D.5b) 


= +me(u - u^) ■ (Ue - 

+ ‘^k{Te - Th){J^^^'> + 
mh 
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